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We evaluate the three bulk viscosity coefficients C11C2 and (,3 in the color-flavor locked (CFL) 
superfluid phase due to phonons and kaons, which are the hghtest modes in that system. We 
, first show that the computation is rather analogous to the computation of the same coefficients 

04 ■ in superfluid ''He, as due to phonons and rotons. For astrophysical applications, we also flnd the 

value of the viscosities when there is a periodic disturbance, and the viscosities also depend on the 
^ , frequency of the disturbance. In a temperature regime that might be of astrophysical relevance, 

we flnd that the contributions of both the phonons and kaons should be considered, and that ^2 is 
much less that the same coefficient in unpaired quark matter. 
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CLc I. INTRODUCTION 

O ' The ground state of three-flavor quark matter at asymptotically large densities and low temperature is the color- 
j-^l flavor locked (CFL) phase Many properties of CFL quark matter have been studied up to now The transport 
or kinetic coefficients of CFL quark matter are of particular interest, as they may allow us to identify whether the 
' CFL phase is realized or not in the core of compact stars, even in the situation where the density is not so large. In 
I> ,. the literature one can find computations of the shear viscosity Q , bulk viscosities H- H , thermal conductivity [3, Q i 
'n]" ' mutual friction @, as well as a general discussion on the kinetics of a CFL superfluid |10|. 

CFL quark matter is a superfluid. The hydro dynamical description of superfluidity differs from that corresponding 
to a normal fluid, as in that case one has to describe the flow corresponding to the quantum condensate, and the 
flow associated to the thermal quasiparticles of the system [HI, Due to the existence of two different flows, more 
\l ■ equations than those corresponding to a normal fluid are needed. In the dissipative regime of the superfluid one can 
also define more transport coefficients than for a normal fluid. In particular, it is known that due to the existence of 
these two different flows one can deflne up to three bulk viscosity coefficients in a superfluid system. 

In this paper we study the three bulk viscosities of CFL quark matter due to superfluid phonons and kaons. These 
collective excitations are the Goldstone modes associated with the spontaneous breaking of baryon and of chiral 
symmetry, respectively. These are the lightest modes in the CFL phase, and thus, one can expect that they provide 
the leading contribution to its transport properties. The contribution to one of the bulk viscosity coefficients due 
■ to kaons only was computed in Ref. Q. The most relevant process for that contribution comes from the decay of 
the neutral kaons into superfluid phonons. In Ref.0] the phonon contribution to the three bulk viscosity coefficients 
was evaluated in a temperature regime where the kaon population was assumed to be thermally suppressed. In that 
case, the most relevant process that contributes to the viscosities is one that changes the phonon number density, and 
implies five phonon collisions. Because the kaon masses are not well-known, it is not clear at which temperatures one 
should consider the kaon contribution into account. 

In this article we compute the three bulk viscosity coefficients of superfluid CFL quark matter due to kaons and 
phonons, assuming a value of the kaon mass gap within a reasonable range of values. While one could naively think 
that the result is the sum of the contributions of the kaons and the phonons, (roughly speaking, adding up the results 
of Ref. Q and Ref.0) we will show that this is not the case, the reason being that kaon decay process also changes the 
phonon number density. We will first show that the computation of the bulk viscosities is totally analogous to that 
corresponding to the phonon and roton contribution to the bulk viscosities in superfluid ^He. Further, and having in 
mind possible astrophysical applications, we will compute the viscosities when there is a periodic disturbance, when 
the viscosities depend also on the frequency of the disturbance. Obtaining the algebraic and numerical values of the 
viscosities is the final goal of this manuscript. 

This paper is structured as follows. In Section |ll] we give a very brief review to the relativistic hydrodynamic 
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equations, and how the bulk viscosities can be interpreted as coefficients that parametrize the deviations of both 
the pressure and chemical potential from its equilibrium values. In Section IIIII we also review the two different 
processes and associated decay rates which contribute the most to the computation of the bulk viscosities. These 
are a five phonon process, and the (electroweak mediated) kaon decay into two supcrfluid phonons. In Section IIVI 
we compute the three viscosity coefficients when there is a periodic disturbance, as a function of the frequency of 
this. In Subsection I V Al we present the analytical results, while the numerical results are displayed in Subsection IV Bl 
Finally, we summarize our main findings in Section [VII We will use natural units h = c ~ kB ~ 1 everywhere in the 
manuscript, except when we show the numerical values of the viscosities. 



II. BULK VISCOSITY COEFFICIENTS IN A RELATIVISTIC SUPERFLUID 



The hydrodynamical equations governing the bulk fluctuations of a non-relativistic superfiuid are essentially different 
from the standard fluid equations. At non- vanishing temperature one has to employ the two-fluid description of Landau 
[ll| , which takes into account the motion of both the superfiuid and of the normal component of the system. In order 
to describe the different dissipative processes one has to introduce more transport coefficients than in a normal fiuid. 
In particular, one has three independent bulk viscosities, as well as the shear viscosity and the thermal conductivity. 
The same occurs in the relativistic version of Landau's two fiuid model, as we briefly review here (see Refs. EBl 
for more explicit details). 

In this paper we will follow the approach derived by Son [ij • This formulation takes into account that superfluidity 
arises after the appearance of a quantum condensate that spontaneously breaks a [/(I) global symmetry. The Gold- 
stone mode associated with this breaking then should appear in any low energy (and low momentum) effective theory. 
Hydrodynamics, being a long spacetimc dynamical formulation, should then necessarily incorporate the Goldstone 
field in their equations. 

The hydrodynamical equations for the superfiuid take the form of conservation laws for both the current, rt^, and 
energy-momentum tensor, T^'^, of the system 

d^^ni; ^ , a^T'"^ = . (1) 

One further adds the Josephson equation, which describes the dynamical evolution of the Goldstone field, or phase 
of the condensate 

uf'd^ip + ^I^o, (2) 

where u^^ is the hydrodynamical velocity, and /i is the chemical potential of the system. The stress-energy tensor and 
the current arc expressed as 

T^'" = (eo + P) u^'u" - Pg'"' + V'^d^ipd'^ip , (3) 
< = n^flu^ - V^d^^ , (4) 

where eg stands for the energy density, P is the pressure, and ^ is a variable proportional to the quantum condensate. 
The energy density obeys the relation 

eo = ST + Hgfl P , (5) 

where S is the entropy of the system. 

Sometimes it is better to write the hydrodynamical equations in terms of the new variable 

w^" = - {d^ip + fiu") , (6) 

which represents a counterfiow momentum, see Ref. In the non-relativistic limit, the spatial component of this 
four vector corresponds to the counterfiow momentum of Landau's two-fluid model formulation [l5| . 

In the absence of dissipation it is possible to show that the entropy current is conserved, i9^(S'm'^) = 0. There is 
generation of entropy, otherwise. The dissipative terms associated to the above hydrodynamical equations have been 
constructed in Ref. [l5|, showing that in the non-relativistic limit they correspond to those appearing in Landau's 
two-fluid model. After defining a comoving frame, where = (1,0,0,0) one can get the terms associated with all 
the viscosity coefficients and thermal conductivity. 

In this article we will only be interested in the three bulk viscosity coefficients, which modify the equilibrium 
hydrodynamical equations. For purposes of interpretation of these coefficients, and for the explicit computations in 
the remaining part of this article, we can see them as follows. The friction forces due to the bulk viscosities can be 
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understood as modifications, with respect to their equihbrium values, in the main driving forces acting on the normal 
and superfluid components. These forces are given by the gradients of the pressure P and the chemical potential fi, 
respectively. One can write 

P = P + ^P = P-Cl^^iV^wn-C2^^u^', (7) 
H = fl + Sfi = fi ~ Csd^^iV^w^ - Cid^uf^ , (8) 

where P and /2 are the equilibrium pressure and chemical potential, respectively. According to the Onsager symmetry 
principle, the coefficients should satisfy Ci = C4i so that in fact there are only three independent coefficients. It is 
also important to stress that the requirement of positive entropy production imposes that <^2 , Ca should be positive 
and that C^f < C2C3 [HI- While (2 has the same meaning as the one that occurs in a normal fluid, Ci and Cs refer to 
dissipation that occurs when the superfluid counterflow is not incompressible. 

The expressions given in Eqs. (|7I8|) are interesting, as they show us that in order to compute the different viscosity 
coefficients we will simply have to evaluate the modifications of both the pressure and chemical potential when the 
system leaves the equilibrium configuration. This is what we will do in the remaining part of this article for the CFL 
phase, using this hydrodynamical formulation, where n^' is the quark current, and fi the quark chemical potential. In 
that case the parameter V'^ can be computed using a Ginzburg-Landau formalism, see Ref. (TB] for the computation 
of V for the CFL superfluid. 



III. PARTICLE NUMBER CHANGING PROCESSES TO RESTORE EQUILIBRIUM AFTER A 

COMPRESSION OR RAREFACTION 

After an expansion (or rarefaction) of a system in equilibrium, the pressure diminishes (or increases above) its 
equilibrium value. Microscopic processes that involve reactions that change the particle densities take place to restore 
the equilibrium configuration. These kind of processes are needed as the system cannot equilibrate to a different value 
of the volume without modifying the density of particles. 

We consider CFL matter at relatively low temperature, where T ^Tc, where Tc is the critical temperature for the 
onset of superfluidity, which is believed to be above several tens of MeV. In this temperature regime we expect that all 
the transport coefficients are dominated by the lightest degrees of freedom of CFL quark matter, as the contribution of 
the heavier modes is Boltzmann suppressed. The lightest modes in the CFL phase are the superfluid phonon, which is 
the Goldstone mode associated with the spontaneous breaking of U{1)b and remains always massless, and the kaons, 
which are pseudo Goldstone modes associated with chiral symmetry breaking. Effective field theories associated to 
these degrees of freedom can be constructed from QCD [l^, Il7l - l20l | . 

In this Section we review the most relevant reaction rates of particle number changing processes that involve the 
CFL light degrees of freedom, and which are relevant for the computation of the bulk viscosity coefficients [1, @ . 

1. Emission, and absorption, of phonons by phonons. The superfluid phonons are massless degrees of 
freedom. Their dispersion law is of the form 

Ep = CsP + Bp^ + -- - , (9) 

where Cs is the speed of sound, and B is also a density dependent parameter that in the CFL phase obeys B < 
[2l| . As realized in Ref. 0, and because B < 0, the first possible phonon number changing reaction involves five 
phonons (f + y:> < — > ip + ip + ip, as the so-called Beliaev process that describes the decay of one phonon in two is only 
kinematically possible for a dispersion law with positive values of B. 

It has been also found that the decay rate associated to the five phonon reaction behaves very differently for large 
or small angle scatterings, being the last processes those which dominate. The decay associated to five phonon small 
scatterings has not been computed, but power counting techniques allows one to estimate it as behaving as 

Tph = aphglgl-^^ , (10) 

where 53 and 54 are the phonon self-coupling constants associated to the three-phonon and four-phonon vertices, 
respectively, and aph is an unknown number. 

Taking into account the value of the strange quark mass, to^, one can find the value of all the phonon self-coupling 
parameters to leading order in a m^Z/i^ expansion. In the density of interest for compact stars, this is the typical 
small parameter used to expand our knowledge of CFL at asymptotic densities. In particular, one finds [H| 

9 1 / ?T!,^ \ n / ml\ -K^ f m'l\ , , 
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Furthermore, a microscopic computation with vanishing quark masses allows one to find the phonon dispersion law 



at cubic order [21 1 



2. Decay of neutral kaons (and antikaons) into two superfluid phonons. This is a process mediated by 
electroweak interactions, which allow for this flavor-changing reaction. First, a kaon is converted into a virtual phonon, 
which finally decays into two other superfluid phonons, so the total reaction can be seen as K'^{K^) i — y (p + (p. The 
neutral kaons/antikaons have a dispersion law of the form 



EkQ/ko = TfJ-KO + V «^P^ + ™K0 ) (13) 

where v = ^, and ii'j^i is an effective chemical potential given by 

eff _ ml -ml 

where rud is the down quark mass, that in this work we will take as negligible in front of the strange quark mass. The 
neutral kaon mass mj^o was computed in the limit of asymptotic high density [isl [l9|. However, at more moderate 
values instantons effects, which are not under full control, give a contribution to all the meson masses [l^, [2^. As a 
result, the kaon mass is poorly known, and we will take it here as an unknown parameter. 

The decay rate associated with this process has been computed in Ref . 3 , where it has been shown that in a good 
range of temperatures it can be approximated by 

^2 f2 f2 ™2 „vp/T 

'^"'^'^ ^ 1^^' + -J^P^" (e-7T _ 1)2 ' (15) 
where the low energy effective field theory coupling constants read 

f2 21~81og(2) fi^ /^v V r ^nwr^ 



where Gi^ is the Fermi constant and 



2 _ e//,2 

t^KO /-|7\ 



KO 

is the momentum where the virtual phonon that mediates the reaction is nearly on-shell. This expression becomes 

r 2 

invalid for very low temperatures, T <C = ^tt- large values of T, which are outside the value of physical 

n 

interest In our numerical computations for the bulk viscosities we will use this approximated form of the kaon 
decay rate, as in the low T regime where this expression is not valid the kaon decay rate becomes meaningless for the 
computation of the frequency dependent bulk viscosities (see Table 2 of Ref. Q)- 

Other processes that involve the kaons are very much suppressed with respect to the one described above. In 
particular, it is possible to show that the rate of the reaction + ip < — > ip is suppressed with respect to the decay 
rate considered above. In Ref [1] , it has been also shown that this is the case for the reaction rates involving charged 
kaons and leptons. 



IV. PHONON AND KAON CONTRIBUTION TO THE BULK VISCOSITY COEFFICIENTS 



In this Section we will compute the phonon and kaon contribution to the bulk viscosities following the same 
strategy devised by Khalatnikov to compute the bulk viscosities coefficients due to phonons and rotons in He*. The 
computation is quite analogous, even in the CFL superfluid the kaons are collective modes with different quantum 
numbers than the phonons, while in He* phonons and rotons are collective modes with the same quantum numbers 
: they are actually the same excitations but seen in a different momentum regime. In any case, all what matters in 
the computation is to have two well-deflned different collective modes which interact with each other. 

The superfluid phonon and the kaon are bosons. At equilibrium the number densities of these particles are entirely 
determined by the temperature T and by their dispersion laws, which depend on the quark chemical potential, and on 
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QCD parameters, such as the gauge coupling constant and quark masses. Small deviations from equilibrium can be 
parametrized with some "fake" phonon and kaon chemical potentials, (5/iph and Sfik, respectively [2^. Equilibrium is 
restored by microscopic processes that modify the particle number densities. The dynamical evolution of the phonon 
and kaon particle densities will be governed by the following equations 

d^^iriphU^) = -iphphSiiph-lphkSiik, (18) 
<9p(7ifeu^) = -JkkS^J.k - IkphSfiph , (19) 

where the different 7^- parameters arc kinetic coefficients which have to be symmetric in their indices, and which can 
be related to the decay rates associated to the microscopic processes under consideration |l2| . If we take into account 
only the scattering processes mentioned in the previous Section one then has 

Tph +4:r K,ph ^^K,ph ^K,ph /„„N 

Iphph = 7^ , Jphk = 7kph = J, , Ikk = rp ■ i^Uj 

As emphasized before, the equilibrium phonon and kaon densities are functions of and T, or cquivalently, 
and S. For small departures from equilibrium, one can follow the change in the phonon and kaon number densities 
by following the dynamical evolution of Uq and S. To this end, we split all thermodynamical and hydrodynamical 
variables into their equilibrium values, denoted with a bar, and a small fluctuation, denoted with a 5. Thus, we write 

Uph = fiph + Suph , (21) 

and similarly for the remaining thermodynamic and hydrodynamical variables. Wc then write 

where the derivatives are taken at the equilibrium values. Because the deviation out of equilibrium is assumed to be 
sufficiently slow we can use the linearized hydrodynamical equations to evaluate the dynamical evolution of Suq and 
6S m 

u^'d^Snq = ^uqd^u^^ ^ d^iV^w^') , (23) 
u^'d^.SS = -Sd^uf" , (24) 

to extract the value of the phonon and kaon chemical potentials. One finds 

S^^p, = 1^ (itd.u'^ - ltd,{V'wn) - -^-^ {I^d.u^ ~ I^d,{V'wn) , (25) 

Ikph ~ Ikklphph V / Ikph ~ Ikklphph 

Sf^k = (itdpu^ - ifd^iV^wn) + {I^d^u'^ - I^d^iV^wn) , (26) 

^kph ^ Ikklphph V ^ Ikph ~ Ikklphph 



where we have defined the functions 

drii _ - drii _ dn 



Il = jr^, r, = n,-S-^-nq-^ , i = ph,k. (27) 



Finally, we compute the change in both the pressure and quark chemical potential caused by the presence of a 
non-vanishing and (5/ip/i 

dP dP 

SP = j: (5Mph + j^Sfik , (28) 

S^i = -p^SfipY, + p^Sfik ■ (29) 

Using the thermodynamical relation Eq. ([5]) and the identity [l^ 

deo = TdS + fiduQ - nphdfiph - Ukd^k , (30) 

one finds 

= t = ph,k. (31) 
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The bulk viscosity coefficients can tfien be obtained by identifying tlie different coefficients in tlie variations of tlie 
pressure and cfiemical potential that are proportional to the divergences of the two velocity vectors, using Eqs (|7l8p . 
One then finds 

Ci = (if + 2/^) (it + 2/|) - -^I^I^ , (32) 

t p/i ^ ' ^ ' t K,ph 

(2 = ^ (if + 21^) (if + 21^) + T^il^r , (33) 

^ph ^ ' ^ ' i K.ph 

(3 ^ ^ (if + 21^) (if + 21^) + -^{II? , (34) 
C4 = Ci ■ (35) 

We note that the last equality is related to the symmetry properties of the kinetic coefRcient jkph ■ If the symmetry 
on the indices of 7^ is violated, then the Onsager symmetry principle is not fullfllled, and then Ci 7^ C4- 

As emphasized at the beginning of this Section, we want to stress that the expressions we found for the viscosities 
are very similar to those found for the phonon and roton contributions to the bulk viscosities to ^He (see Ref. [HI). 

In the situation where the temperature is much smaller than the kaon gap T <C Sm = — yU^o, then the 
contribution of the kaons is thermally suppressed. In Ref. @ the three bulk viscosity coefficients due to the phonons 
were evaluated, that is, in the situation where I^ = 0. Because the kaon gaps are poorly known, it is not clear at 
which values of T the kaons might become relevant or not for transport in the CFL phase. 

While the phonons are never thermally suppressed because they are massless degrees of freedom, one can still find 
a regime where the kaons dominate the contribution to the bulk viscosities. This may happen when Tph 3> ^K,phi as 
then the first terms of Eqs. p2ll34p would be suppressed. Because of the high power dependence of the phonon decay 
rate, one may expect that this is the situation that one may encounter at high temperatures. 

One important lesson deduced from the expressions written above is that in the domain where the kaons are not 
thermally suppressed the pure phonon and kaon contributions cannot be disentangled, as the bulk viscosities contain 
products of both phonon and kaon thcrmodynamical functions. That is, the total bulk viscosity in the presence of 
phonons and kaons is not simply the sum of the phonon and kaon contributions. 



V. PHONON AND KAON CONTRIBUTIONS TO THE FREQUENCY-DEPENDENT BULK 

VISCOSITIES 

For astrophysical applications it is important to find the value of the bulk viscosities when there is perturbation 
which is periodic in time, with an angular frequency coc- Then, one also assumes that the time evolution of all the 
physical variables is also periodic. The bulk viscosities then turn out to be also a function of Wc, and in this Section we 
will compute the frequency-dependent value of these transport coefficients. In Subsec. FV Al we present the algebraic 
expressions of the viscosities, and we leave the presentation of the numerical results for Subsec. IVBI 



A. Analytical results 

The evaluation of the frequency dependent bulk viscosities when only phonons arc considered was done in Ref. 0] . 
In Ref. [3, C,2{'^c) due to kaons was also computed. Here we want to get the values of the frequency dependent 
viscosities when both the phonons and the kaons are taken into account. It is not a priori obvious how the coefficients 
will behave, as we have already seen that in a non-periodic situation one cannot simply sum the phonon and kaon 
contributions to the bulk viscosities. 

We will work in the comoving frame, where — (1,0,0,0). Furthermore, as in Ref. [3| we will compute the 
viscosities assuming that the temperature is constant. In this regard, our computation here differs from that carried 
out for the phonons in Ref. @ , where such an assumption was not considered. 

For purposes of comparison with the results of a non-periodic case presented in the previous Section, and the results 
obtained in Ref. we will first find the viscosities neglecting the variation of the phonon number, that is, we will 
assume momentarily that Suph = S^ph = 0. We thus use only the kaon and quark density equations. In the comoving 
frame, and taking into account the explicit time dependence of the variables, one has 



iuJc Suk +nk divu = ~jkk Sfik , 
iujc Suq +nq divu + div(y^w) = 



(36) 
(37) 
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Wc now express the fluctuations of the particle densities in terms of a change in both the quark and kaon chemical 
potentials 

5n, = ^-^5^,+ ^5„u. (38) 
5nk = Sfi + Suk , (39) 

where the derivatives of the particle number densities are evaluated in equilibrium, where one has = 0, see 
Appendix |X]for explicit details. In this way, we can solve for 5^ and 5^k as a function of the different thermodynamical 
variables and the angular frequency. The solutions, and therefore also the bulk viscosities, are complex quantities. 
However, for the computation of the dissipation in energy only the real part of these coefficients is needed, so we will 
only extract the real part of the viscosities. 

Using Eqs. (|7l8p and Eqs. p8l29|) . we find that the real part of the viscosities can be expressed as 

C.(wc) - , (40) 



where 



1 1 /_ dn„ _ duk 



1 1 (dnk\ 



2 



2 



dp. 



and 



and we have defined the variable 



_ duq duk duq duk 



d^ dfik d^k 



(45) 



It is interesting to note that the condition for positive entropy production is saturated, and we have (1 = C2C3- 
With the results presented above, we can check that the formal expression of C2 agrees with the bulk viscosity 
computed in Ref. Q , as one has = (see Appendix . We can also check that in the limit we obtain 

the values of the bulk viscosities obtained in the previous Section, in the situation where one assumes no thermal or 
entropy variation, and for /f'* = J^''' = 0. In order to check this limit, one simply has to realize that at constant 
temperature 

to write down all the pieces Ci / jef / in terms of the decay rate and also the and /| functions introduced in Eq. (|27p . 

Let us consider now the presence of fluctuations in the phonon sector. The set of linear equations that wc have to 
solve should now also involve the phonon number density equation. The system of equations to solve is 

iujc Suph +nph divu = - jphph Sfiph - Iphk Sfik , (47) 
iuJc Suk +nk divu = -jkk S^J,k - Ikph Stipt , (48) 
divu + div{V^w) = . (49) 

One can find solutions for 6fj,, 6fj,k and S^ph if we express the fluctuations of the particle densities in terms of the 
fluctuations of the chemical potentials 

dUq dUq dUq 

5n, = (50) 
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Suk = -7:—S^J.+ ^ — dfik , 51) 

Suph = -TT^Sfl + -;r-^S^ph ■ (52) 

dfiph 

Then, using Eqs. (|7l8p and Eqs. (|28l29p . one can find the real part of the viscosities. In this case we find that the 
frequency dependence of the coefficients is different from the case where only the kaons (or phonons) arc considered. 
The real part of the three coefficients can now be written as 

C.(-c) = , (53) 



where 



O+Ruj^ 



°^T(ifj ■ f^"' 

R = ^ (-^L + ( 1^1 , (55) 




and Di are polynomials of second order in lJc, expressed as 



(56) 



A(wc) = {h O + ujI {b, R + a, Q)) . (57) 

The coefhcients Ci, and bi are different for each bulk viscosity, and depend on both the phonon and kaon thermody- 
namical variables, and on the the decay rates. More explicitly, the coefficients for the first bulk viscosity are 



duph dru 



( n ^Ih. _ n \ _ jph driph duk dug 



duph dfj. 



jph ^K,ph duq f duph dnu \ 



The coefficients for the second and third bulk viscosities are 

V 9^ dii ) l"^^ dji djjL V 9/i d^iph dfi d^ph ) 9^ d^iph ' 

, TPh ^"g f - r , - ^^ph ^"g - 9^ph duk 
+ I2 -Tr~ \nphL + Uk — t; rig — — - — 

, ^K,ph / , _ ^ 9?!, _ f duk duph , , , . 

= -^[i'^rik+npn)^-ng[2— + ^]) (61) 



duq _ dnk\ j-phdng\ Tph (_ duq _ duk 



(62) 



and 



dn„h \ ^ dnk ( duk \ ^ 9n„ 



^ \ dpL J d^k V di-L J dfiph ' 

^Kph f . drvk f duph drik\ duph dug \ 1^ / dnk\ ^ , , 

T \ d^l \ d^^ ^ dti J ^ dfi d^lph) t [dti J ' ^ ' 
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respectively. Further, we have found Ci = C4- 

We have checked that our expressions in the hmit ~^ reproduce the results of the bulk viscosities in the limit 
when there arc no thermal or entropy fluctuations, obtained in the previous Section. While it is not easy to verify 
whether the expressions for the three bulk viscosities are consistent with the constraint of positive entropy production, 
we have checked numerically that this is the case, as one always gets Ci < CaCs- 

Because of the rather non-trivial dependence of the expressions of every bulk viscosity coefficient on the phonon 
and kaon thcrmodynamical variables it is not easy to infer how the viscosities scale with the temperature, or with 
other relevant scales of the problem, such as /i. It is also not easy to sec when the phonons or the kaons give a 
subleading contribution. A numerical analysis is then mandatory. We present the numerical results of the viscosities 
in the following Subsection. 



B. Numerical results 




FIG. 1: Temperature dependence of the first and third bulk viscosities for the CFL quark matter phase due to tliermal phonons 
and kaons for a frequency uic/2ti — 1 Khz. The phonon decay factor is fixed at aph = 1, and we plot here the value of the 
viscosities for five different values of the kaon mass gap. 

We have numerically computed the values of the three bulk viscosities, and we present the results in different 
logarithmic plots for fixed values of the frequency, and as a function of the temperature. We consider temperatures 
in the range of 0.01 MeV <T < 10 MeV. Lower temperatures might be relevant in astrophysical scenarios. However, 
it was found out in Ref. 0] that for lower T the phonon mean free path is much larger than the radius of a compact 
star, and a hydrodynamical description of the phonon fluid ceases to be valid. 

In our numerical computations we have chosen the value of the quark chemical as /i = 400 MeV, the quark gap 
A = 25 MeV and the strange quark mass nig = 120 MeV. We leave as free parameters the kaon gap Sm = mj^o — fi'j^i , 
and aph, the numerical factor present in the phonon decay rate, see Eq. (jlOp. Notice that while we worked using 
natural units in the main part of the manuscript, we use SI units to show the numerical values of the viscosities. It 
is important to realize that the three bulk viscosities have different scalings in length (see the Appendix [B|) . 

The phonon and kaon thcrmodynamical variables that are needed for the viscosities can be computed from the free 
energy of the system (see Appendix \X\ for details), and we use the decay rates shown in Section Hill 

In Fig. [Hand Fig. [2] left we present the numerical values of CiiCs ^-nd (2, respectively, for a fixed value of the 
frequency, lUc/2tt = IKhz, and for five different values of the kaon mass gap, ranging from 6m = 0.1 MeV to 10 MeV. 
We have fixed in these plots the value of aph = 1. Roughly speaking, the plots for the three viscosities present a 
similar shape, that is, a similar dependence on T, although the numerical values for each viscosity are very different 
due to the different scaling dimensions of each coefficient. 

In Fig. [2] right we also present the value of C2 when only the kaon contribution is taken into account, that is, with 
the results of Ref. 3 , or when only the phonon contribution is taken into account (the results here differ from those 
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FIG. 2: Left: Temperature dependence of second bulk viscosity for the CFL quark matter phase due to thermal phonons and 
kaons for a frequency cjc/Stt — 1 Khz, with the phonon decay factor is fixed at aph ~ 1, and for five five different values of the 
kaon mass gap. Right: The same plot but for ^2 due to kaons only, or to phonons only. 



of Ref. Q because here we are doing the computation assuming that the temperature remains constant). The same 
plots for the kaon contribution only for ^1 and Cs present a rather similar shape, so we won't show them here. 

In Fig. 12] we show the sensitivity of one of the bulk viscosities to the value of Uph , the constant in the phonon 
decay rate Eq. (fTO)) and for two fixed values of the kaon mass gap. We can observe that the numerical values of the 
viscosities are rather insensitive to the values of ranging form 10~^ to 10^, except in the high T limit. This is due 
to the high power dependence on T of the phonon decay rate, Eq. (|10p. For this reason, we have fixed the value of 
o-ph = 1 in most of the plots we are showing. 

In Fig. m we show the value of C2 varying the value of the frequency, for Uc = 1, 1/3, 0.1 Khz, and we have shown 
also in the same plot the same results for the kaon contribution only. 

Some generic comments can be formulated from our numerical studies. At very small temperatures the phonons 
dominate the contribution to the bulk viscosities, as the kaons are then thermally suppressed. This fact is not always 
shown clearly in the plots as this tends to occur at T ^ 0.01 MeV, which is outside the temperature regime where 
one may expect the phonons to be in a hydrodynamical regime inside the star. At very large temperature T > 10 
MeV the kaon contribution seems to dominate the behavior of the viscosities. The reason comes from noting that in 
that regime Tph becomes very large, and suppresses the phonon contribution (see comments at the end of Sec. lIV|) . In 
the intermediate temperature regime we always find an interval where first the kaons dominate, followed by a generic 
two-peak/hump structure, where both phonon and kaon contributions are relevant. The temperatures where this 
happens depends on the value of the kaon mass gap. 

One of the main conclusions derived from our numerical studies is that in the regime where the hydrodynamical 
description might be valid for the superfluid phonons, which should be in the range T > 0.01 MeV, and for a range of 
values of kaon gaps that do not exceed the few MeV, both kaons and phonons contribute to the value of the viscosities. 
The superfluid phonon would only be highly dominant for T <^ 0.01 MeV, but as already mentioned, in this regime 
in a compact star the phonons do not collide enough times to maintain a fluid description in that case. Then the 
phonons rather behave as a gas, and one should consider other dissipative effects in the system, such as that coming 
from the collisions of phonons with superfluid vortices 0] . 

In Ref. [1] it was realized that the value of C2('^c) as due to kaons only is much less than the value of the same 
coefficient for unpaired quark matter. From our numerical results, we can say that if we include both the phonon and 
kaon contribution (see Fig. [2] and FigU]) the reduction is still more severe, although this depends on the actual value 
of (5m and on the values of T. 
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FIG. 3: Temperature dependence of the third bulk viscosity, for a fixed frequency lJc = 1 Khz, for two different values of the 
kaon mass gap, and varying the value of the phonon decay factor from Uph = 




FIG. 4: Temperature dependence of the second bulk viscosity for three different values of the frequency, taking Uph = 1, and 
for two different values of the kaon mass gap. For comparison we have also plotte in the same graphs the values of Cii'^c) when 
only the kaon contribution is taken into account. 



VI. OUTLOOK 



We have presented a computation of the frequency dependent bulk viscosities in the low temperature regime of CFL 
quark matter, extending previous results already obtained in the literature. The main contributions to the viscosities 
come from both kaons and superfluid phonons. The main results of this manuscript are presented in Eqs. (|53ll64p and 
Figs. ID to H 

One of the main conclusions of our study is that in the regime where the hydrodynamical description might be valid 
for the superfluid phonons, which should be in the range T > 0.01 McV, and for a range of values of kaon gaps do 
not exceed the few MeV, both kaons and phonons contribute to the value of the viscosities, and one cannot discard 
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any of the two quasiparticle contributions. Our results also agree with the main finding of Ref. [j], that is, that the 
coefficient ^2 in the CFL phase is several order of magnitude smaller than the same coefficient in unpaired quark 
matter. 

The studies of transport coefficients are important in order to find possible signatures of quark matter in astro- 
physical scenarios. The fate of different set of oscillations modes of compact stars are governed by the value of the 
transport coefficients in the star. In particular, the fate of the r-modes of the star [24|] depend on the dissipative 
effects in the star. Previous analysis of the damping of the r-modes have been carried out in the literature [ol. [25l - [27| . 
The results of our study should be taken into account in a more refined study of the the damping of r-modes, where 
the value of the three bulk viscosity coefficients should be relevant 
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Appendix A: Phonon and kaon thermodynamical variables 

The phonon and kaon thermodynamical variables needed for the computation of the bulk viscosities can be extracted 
from the free energy of the CFL system, which reads 

Sfik, Sfiph, T) = T) + Vtkin-, S^k, T) + ilphip-, Sfiph,T) , (Al) 

where fig is the quark contribution, while fl^ and flph are the kaon and phonon contributions, respectively. The 
free energy of the phonons and kaons in the presence of the "fake" chemical potentials that characterize the out of 
equilibrium state is 



rp />OC 

Ui{fi,Sm,T) = J dpp'^ \n {1 - cxp{- {Ei - Sfj.i)/T)) , i = k,ph 



(A2) 



The /i-dependence in the contribution to the free energy of the phonons and kaons appears because the energy of 
the quasiparticle depends on fi. In particular, for the phonons both the speed of sound Cs and the parameter B that 
enters in their dispersion relation, see Eq. (|12p . depend on fi. For the kaons, both their velocity v, mass mxo, and 
chemical potential f/j^o depend on fi, although here we will only consider the dependence of ^^o^, as the dependence 
of the other variables is unknown. 

The different particle densities are obtained from the free energy of the system 

ni = -——, i^q,k,ph. (A3) 

From the free energy of the system we can also compute all the derivatives that appear in the expressions of the 
bulk viscosities. These are computed from Ea. (|Aip . putting S^i = at the end of the computation. 

The phonon thermodynamical variables can be computed analytically. For the computation of the bulk viscosities 
that was done in Ref. [61 it was realized that the phonon dispersion law at non-linear order was needed. The explicit 
expressions of the thermodynamical functions needed in our computation can be found in that reference. The kaon 
thermodynamical variables have to be computed numerically. 

It is interesting to note that the following conditions are satisfied 

dnk _ dnq_ duph _ dug duph _ duk _ ^ ^^^^ 

d^k ' dfj. dfj.ph ' dfj.k dfj,ph 

Further, we will use the relation m 



dnk _ ml duk 
djjL 2/i2 d^k 



(A5) 
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Appendix B: Dimensions of the bulk viscosity coefficients 

The bulk viscosity coefficients introduced by Gusakov in Son's formulation of the relativistic superfluid hydrody- 
namics arc related to those introduced by Khalatnikov by mass factors. More explicitly, one has 

/-Kh ^Kh A /-Kh Cs ,rn\ 

m 771^ 

If M, L, T refer to scales of mass, length and time, respectively, then the three bulk viscosity coefficients in Gusakov's 
theory have the following dimensions 

[Cl] = ML^T-^ , [C2] = ML-^T-^ , [C3] = ML^r-i . (B2) 

In all the main part of the paper we worked using natural units. The plots are given in units of SI. For the conversion 
of the values of the bulk viscosities from natural units to IS units it is useful to remember that 

lMeV= 1.78- 10-2^g , IMeV = 1.97- 10~"cm , 1 MeV = 6.58 • lO'^^sec . (B3) 

When the bulk coefficients are expressed in natural units they seem to behave almost similarly, that is, one can see 
that they have a similar T dependence. In the SI system, and due to the different scalings in length of the coefficients, 
the coefficients show a numerical behavior rather different, as shown in the plots of Fig. [l]and Fig. [51 the ultimate 
reason being the different scaling in L of every coefficient. 

We note here that in the analysis of the r- modes of the GFL phase done in Ref. ^27|], it was claimed that only Ca 
is important for the computation of the damping of the r-modcs. However, the dimensions of Ca in that manuscript 
have been wrongly assumed to be the same as those of C2, which probably affects the final numerical results of the 
whole analysis when converting the different magnitudes from natural units to SI units. 
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